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Abstract 



A product of gauge fields generated by the Yang-Mills gradient flow for 
ed ■ positive flow times does not exhibit the coincidence-point singularity and a 

local product is thus independent of the regularization. Such a local product 
Q . can furthermore be expanded by renormalized local operators at zero flow 

r^, ! time with finite coefficients which are governed by renormalization group 

equations. Using these facts, we derive a formula that relates the small flow- 
time behavior of certain gauge invariant local products and the correctly- 



cn ' normalized conserved energy-momentum tensor in the Yang-Mills theory. 

Our formula provides a possible method to compute correlation functions of 
O ! a well-deflned energy-momentum tensor by using the lattice regularization 

Tjj- \ and the Monte Carlo simulation. 

Keywords: Energy-momentum tensor, Yang-Mills gradient flow. Lattice 
gauge theory, Wilson flow 



1. Introduction 

Although the lattice regularization provides a very powerful non-perturbative 
formulation of fleld theories, it is unfortunately incompatible with fundamen- 
tal global symmetries quite often. The most well-known example is the chi- 
ral symmetry [l|, |2| and the supersymmetry is another infamous example |3| 
needless to say the translational invariance. When a regularization is not 
invariant under a symmetry, it is not straightforward to construct the cor- 
responding Noether current that is conserved and generates the symmetry 
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transformation through Ward-Takahashi (WT) relations. This makes the 
measurement of physical quantities related to the Noether current in a solid 
basis very difficult. To this problem, one can imagine at least three possible 
approaches. 

The first approach is an ideal one: One finds a lattice formulation that 
realizes (a lattice-modified form of) the desired symmetry. If such a for- 
mulation comes to hand, the corresponding Noether current can easily be 
obtained by the standard Noether method. The best successful example in 
this sort is the lattice chiral symmetry j^, |5|, la, lD, Isl, M, UM, luj that can be 



defined with a lattice Dirac operator that satisfies the Ginsparg- Wilson rela- 
tion [12]. Although this is certainly an ideal approach, it appears that such 
an ideal formulation does not always come to hand, especially for spacetime 



symmetries (see Ref. [13|, for example, for a no-go theorem for the super- 
symmetry). 

The second approach is to construct the Noether current by tuning coef- 
ficients in the linear combination of operators that can mix with the Noether 
current under lattice symmetries|j For example, for the energy-momentum 
tensor — the Noether current associated with the translational invariance and 



rotational and conformal symmetries [ij, ll5|| — one can construct a conserved 



lattice energy-momentum tensor by adjusting coefficients in the linear com- 
bination of dimension 4 operators 16|, llTI j; the overall normalization of the 



energy-momentum tensor has to be fixed in some other wayo Although 
this method is in principle sufficient when the energy-momentum tensor is 
in "isolation", i.e., when the energy-momentum tensor is separated from 
other composite operators, as in the on-shell matrix elements, it is not 
obvious a priori whether one can control the ambiguity of possible higher 
dimensional operators which may contribute when the energy-momentum 
tensor coincides with other composite operators in position space. This 
implies that it is not obvious whether the energy-momentum tensor con- 
structed in the above method generates correctly-normalized translations 
(and rotational and conformal transformations) on operators through WT 



-'^Here, we assume that the fine tuning of bare parameters to the target (symmetric) 
theory is made. 

^A somewhat different approach on the basis of the A/" = 1 supersymmetry has been 



given in Refs. J18lll9|. 

•^It might be possible to employ the "current algebra" for this, like for the axial cur- 
rent El. 



relations. (If the energy-momentuni tensor generates correctly-normalized 



translations, it is ensured UM (see also Sec. 7.3 of Ref. [22|) that the trace 
or conformal anomaly |23|, |2J] is proportional to the renormalization group 
functions B H HO 

The third possible approach is to utilize some ultraviolet (UV) finite quan- 
tity. Since such a quantity must be independent of the regularization adopted 
(in the limit in which the regulator is removed), there emerges a possibility 
that one can relate the lattice regularization and some other regularization 
that preserves the desired symmetry. This methodology can be found for 
example in Ref. |28| (see also Ref. [29|), where an ultraviolet finite represen- 
tation of the topological susceptibility is derived. Although the derivation of 
the representation itself relies on a lattice regularization that preserves the 
chiral symmetry j^, Isl Ig, |7|, Isl. |9| . llOl . |ll|, one can use any regularization (for 
example, the Wilson fermion |30|) to compute the representation because it 
must be independent of the regularization. 

In the present paper, we consider the above third approach for the energy- 
momentum tensor, by taking the pure Yang-Mills theory as an example. For 
this, we utilize the so-called Yang-Mills gradient flow (or the Wilson flow in 
the context of lattice gauge theory) whose usefulness in latti ce g auge theory 



has recently been revealed [31|, l32|, |33|, |3J, |35|, |36|, l37|, l38|, |39 



feature of the Yang-Mills gradient flow is its robust UV finiteness [33 



A salient 
More 



precisely, any product of gauge fields generated by the gradient flow for a 
positive flow time t is UV finite under the standard renormalization. Such a 
product moreover does not exhibit any singularities even if some positions of 
gauge fields coincide. The basic mechanism for this UV finiteness is that the 
flow equation is a sort of the diffusion equation and the evolution operator 
in the momentum space ~ e"*'^ acts as an UV regulator for t > 0. This 
property of the gradient flow implies that the definition of a local product 
of gauge fields for positive flow times is independent of the regularization. 
In our present context, there is a hope to relate quantities obtained by the 
lattice regularization and obtained by the dimensional regularization with 
which the translational invariance is manifest. 
As noted in Ref. 



33 



on the other hand, a local product of gauge fields 
for a positive flow time can be expanded by renormalized local operators of 
the original gauge theory with finite coefficients. Those coefficients satisfy 
certain renormalization group equations that, combined with the dimensional 
analysis, provide the information on the coefficients as the function of the flow 
time. Because of the asymptotic freedom, one can then use the perturbation 



theory to find the asymptotic behavior of the coefficients for small flow times. 
By using above properties of the gradient flow, one can obtain a formula 
that relates the small flow-time behavior of certain gauge invariant local 
products and the energy-momentum tensor defined by the dimensional regu- 
larization. Since the former can be computed by using the Wilson flow with 



the lattice regularization [SfJ, l32|, |33|, |3J, |35|, |36|, |37|, |38|, |39[ and the latter 
is conserved and generates correctly-normalized translations on composite 
operators, our formula provides a possible method to compute correlation 
functions of a correctly-normalized conserved energy-momentum tensor by 
using the Monte Carlo simulation. 

In the present paper, 
unless otherwise stated. 



we follow the notational convention of Ref. 33 



2. Yang— Mills theory and the energy-momentum tensor 

2.1. The energy-momentum tensor with the dimensional regularization 

In the present paper, we consider the SU{N) Yang-Mills theory defined 
in a D dimensional Euclidean space. The action is given by 



from the Yang-Mills field strength 



S=^,jd^xF;,{x)F;^{xl 



(2.1) 



F^,u{x) = d^A^ix) - d^A^ix) + [A^{x),A^{x)]. 



(2.2) 



We set 

D = A-2e, (2.3) 

and then the mass dimension of the bare gauge coupling go is e. 

Assuming that the theory is regularized by the dimensional regulariza- 
tion (for a very nice exposition, see Ref. J40|), one can define the energy- 
momentum tensor for the system (12. ip simply by (see, for example, Ref. |4l|) 



-^ fj,u\X ) 



9l 



F;,{x)F^.,{x) - -5,.F^A^)F;^ix) 



(2.4) 



up to terms attributed to the gauge fixing and the Faddeev-Popov ghost 
fields, which are irrelevant in correlation functions of gauge invariant opera- 
tors. Note that the mass dimension of the energy- momentum tensor is D. 



The advantage of the dimensional regular izat ion is its translational invari- 
ance. Because of this property, the energy-momentum tensor naively con- 
structed from bare quantities, Eq. (12 ■4p . is conserved and generates correctly- 
normalized translations through a WT relation, 

Jd^'x {d,T^,{x)0) = - {d,0) , (2.5) 

where it is understood that the derivative in the right-hand side is acting 
all positions in a gauge invariant operator O. Used in combination with the 
dimensional counting and the gauge invariance, this WT relation implies that 
the energy- momentum tensor T^y{x) is finite 42|, |26| and thus in the minimal 
subtraction (MS) schemep 

T,.(x)-(r^.(x)) = {T,4^(x). (2.6) 

The finiteness of the energy-momentum tensor (12.41) provides further use- 
ful information on the renormalization of dimension 4 gauge invariant opera- 
tors. The gauge coupling renormalisation with the dimensional regularization 
is defined by 

gl ^ fi'^g'Z, (2.7) 

where /i is the renormalization scale and Z is the renormalization factor. In 
the MS scheme, 



1 

1 -- 
e 



bo9' + lhg' + 0{g') 



+ 0(1), (2.8) 



and 

'» =48^- '■ = 384^- P-8' 

From the rotational invariance that the dimensional regularization keeps, we 
see that the operator-renormalization possesses the following structures|j 

= Zt {F^,F:^}^ (x) + Zm6,. {F^.F^^}^ (x), (2.10) 



Here, we define the renormalized operator by subtracting its vacuum expectation 
value. In the perturbation theory using the dimensional regularization, this subtraction is 
automatic. 

^Here again, we define renormalized operators by subtracting their vacuum expectation 
values. 



and 



f^Mf;a^) - {f;a^)f^,A^)) = zs {f;^f;^]^ (x). (2.11) 

Substituting the above relations in Eqs. (12 ■4p and (12.61) . we have 



1 

72 



^/^-"^-^ 
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Zt {f::,f^%}^ (x) -^{Zs- ^Zm)5,. {f;.f;^]^ 



\x 



(2.12) 



Since the left-hand side is finite for e — > 0, in the MS scheme in which 
only pole terms are subtracted, we infer (by considering the cases, \i i^ v 
and /i = z/), 

Zt = Z = 1-6o/- + 0(/), (2.13) 



and 



4Z 



M 



z. 



(2.14) 



2.2. Implications oj the trace anomaly 

Another important property of the energy- momentum tensor (12. 4p is the 
trace anomaly |25|, |26|, |27| , 



d^u i-'^ ^^.ufJl [X) 



/3 



2r 



{F^.F^.}^{x). 



(2.15) 



By Eq. (12. 6p . this relation is equivalent to 

S,. [T,^{x) - (T,.(x))] = e-^F;^(x)F;,(x) - (e^F;:,{x)F;^{x) 



e^O 



P 



{F^.F^.}^ix). 



2n3 L P'^ P^JR 

In Eqs. (I2.15P and (I2.16p . /3 denotes the /3 function for D = 4, defined by 



(2.16) 



/3 



d_ 



'^'^^l^)o''^' 



(2.17) 



where the subscript implies that the derivative is taken while the bare 
quantities are kept fixed. Eqs. (12. 8p and (12. 7p yield 



(3 



-bog'-hg' + 0{g' 



(2.18) 



Then substituting Eqs. (ETT^ into Eq. ^HM and using Eqs. (ITT^ and (I^THjl . 
we observe that 



5pA {i^;.i^t }^ (x) = (i - ^) {C^;.} 



i? 



(x), (2.19) 



that is, the contraction with the metric and the minimal subtraction, the 
subtraction of 1/e poles, do not commute; this is a pecuhar but legitimate 
property of the dimensional regularization 40| . 



Also, substituting Eqs. (ET]) and fl2AT|) into Eq. ( 1216|) . we see 

6% ^ -^. (2.20) 

z g 

In the MS scheme in which only pole terms are subtracted, this implies 

t^A r. . , ^. 4^ 



Zs= (l-^-)Z=l + 0(/), (2.21) 



and Eq. (12.141) then shows 



/5 ^1 ^0 2^ , D^ 4.\ 



Zm = -^Z- = -fg'- + 0{g'). (2.22) 

4(7 e 4 e 

We thus observed that all the renormalization constants in Eqs. (I2.10p and (12.111) . 
Zt, Zm and Zs, in the MS scheme can eventually be expressed by the gauge 
coupling renormalization constant Z in Eq. (12. 7p . 

3. Yang— Mills gradient flow and the small flow-time expansion 

The Yang-Mills gradient flow deflnes a -D + 1 dimensional gauge poten- 
tial B{t,x) along a flctions time t, according to the flow equation, 

dtBf,{t, x) = D^G^^{t, x) + aoDf,d^B^{t, x), (3.1) 

where the D + 1 dimensional field strength and the covariant derivative are 
defined by 

G^,(t,x) = d^B,it,x)-d,B^{t,x) + [B^{t,x),B,{t,x)], (3.2) 

and 

D^ = d^ + [B^,l (3.3) 



respectively. The initial condition for the flow is given by the D dimensional 
gauge potential in the previous section: 

B^{t = Q,x) = A^{x). (3.4) 

In Eq. (13.11) . the last term is introduced to suppress the evolution of the field 
along the direction of gauge degrees of freedom. Although this term breaks 
the gauge symmetry, it does not affect the evolution of any gauge invariant 



operators [31] ■ Note that the mass dimension of the flow time t is —2. 

Now, from the field strength extended to the -D + 1 dimension (13. 2p . we 
define a. D -\-\ dimensional analogue of the energy-momentum tensor by 

f/^,(t,x) ^ Gl^{t,x)Gl^{t,x) - \5,,Gl,{t,x)G^^,{t,x). (3.5) 

Although this is similar in form to the original energy-momentum tensor (12. 4p . 
it is not obvious a priori how this -D + 1 dimensional object and Eq. (12. 4p are 
related (or not). To find the relationship between them is the principal task 



of the present paper. We also use the density operator studied in Ref. [31 



E{t,x)^\G%{t,x)G%{t,x). (3.6) 



Now, as shown in Ref. [33|], for t > 0, any correlation function of B^{t, x) 
is UV finite after the standard renormalization in the D dimensional Yang- 
Mills theory. This property holds even for any local products of B^{t,x) 
such as Eqs. ( 13. 5p and ( 13. 6p . Also, for small flow times, a local product 
of -B^(t, x) can be regarded as a local field in the D dimensional sense because 
the flow equation (13. ip is basically the diffusion equation along the time t 
and the diffusion length in x is \/8t. These properties allow us to express. 



as explained in Sec. 8 of Ref. [33|, U^u{t,x) and E{t,x) as an asymptotic 
series of D dimensional renormalized local operators with finite coefficients. 
Considering the gauge invariance and the index structure, for D = 4, we can 
write 

f/^.(t,x) =CT(t){T^.}^(x) + C5(t)v{^^;.^;.} {^) + 0{t), (3.7) 

where abbreviated terms are the contribution of operators with the mass 
dimension higher than or equal to 6. For Eq. ( 13. 6p . we similarly have 

E{t,x) = {E{t,x)) + CE{t)^^F;,F;}^ {x) + 0{t). (3.8) 



csit) = TTlMt), (3.10) 



We note that, when the renormahzed gauge couphng is fixed, U^^it^x) fl3.5p 
is traceless for D = 4, 

5^.f/^.(t, x) = 2eE{t, x) ^ 0, (3.9) 

because E{t,x) f l3.6p is finite [31i] and does not produces a 1/e singularity 
(this explains why there is no c number expectation value term in Eq. (13. 7p ). 
Thus, considering the trace part of Eq. (13. 7p . we see that the coefficients 
CT(t) and cs(t) are not independent and related by, for D = 4, 

A 

because of the trace anomaly (I2.15P . 

By eliminating the renormahzed action density from Eqs. (13. 7p and (13. Sp . 
we have 

{T^J^(x) = -^U,,{t,x) - -^0L_6,, [E{t,x) - {Eit,x))] + 0{t). 

(3.11) 
This expression relates the energy-momentum tensor (12.61) and the short fiow- 
time behavior of gauge-invariant local products defined by the gradient fiow. 
Thus, once the coefficients are known, one can extract the energy-momentum 
tensor from the t — ?■ behavior of the combination in the right-hand side. 

4. Renormalization group equation and the asymptotic formula 

4-.1. Renormalization group equation for the coefficients 
We now operate 

on the both sides of Eq. (13. 7p . Since the left-hand side of Eq. (13. 7p . i.e., 
Eq. (13. 5p . is entirely expressed by bare quantities through the fiow equa- 
tion (13. ip and the initial condition (13. 4p . the action of (14. ip on the left-hand 
side identically vanishes. In the right-hand side, this vanishing-ness must 
hold in each powers of t. Thus we infer that 

{^^) ^T{t){T,,}^{x) = 0, (4.2) 



i)Ml- 



-"— ) csit)\-F;^F;^\ (x)=0. (4.3) 



For the first relation f l4.2p . we recall that the energy-momentum tensor 
is not renormalized as Eq. fl2.6p . Then by expressing the operation (14.1 p in 
terms of renormalized quantities, we have 



For Eq. (14. 3p . on the other hand, from Eq. (12. lip . 



'*l+4^''^'''^'*>=°' 



where 



75 






Eqs. (E2ID, (EZD and (IZTTj) yield 



d f 13 



73 = -9 , . 3 
Similarly, for Eq. ( 13. Sp . we have 



d 8 

f^TT + P 



d 



dg 
d 



2M' + 0(/ 



(E(t,x)) = 0, 



^a;i+^a^+^^'^-(^) = °' 



and thus 



5/i (95-7 CE{t) 



0. 



(4.4) 

(4.5) 
(4.6) 

(4.7) 

(4.8) 
(4.9) 

(4.10) 



By the standard argument and from the fact that dimensionless quantities 
can depend on the renormalization scale /i only through the dimensionless 
combination y/Stfj., the above renormalization group equations imply that 



CT{t){g; fj^) = CT(to)(^(-0;/^o 
cs{t){g; fJ') = exp 



dClsim')) 



Cs{to){g{-^);fio), 



e {E{t,x)) (g-fi) = tl {E{to,x)) (^(-0;/io) 

cs{t) Csjto) 

-\9\N = — 7T\\9{-0\M^ 



ceH) 



CEita 



(4.11) 
(4.12) 
(4.13) 
(4.14) 



10 



where the dependence on the renormahzed gauge couphng and on the renor- 
mahzation scale has been exphcitly written. In these expressions, the running 
couphng g{^) is defined by 

^ = /3(M0), m) = 9, (4.15) 

and we introduced a variable 



^ = ln^B^. (4.16) 



In the one-loop order, the running couping (14.151) is given by 

^^"^^ ~ 26;; -e + 1/(2^2) - Wo^h^(^Mf^HHV^^' '^ 



where A is the A parameter in the one-loop level, 

A = ^e-i/(2^o3'), (4.18) 

and the integral appearing in Eqs. (14.121) is 

/ ^ de is m')) = ^1^ [g' - 9i-0'] ■ (4.19) 

Jo "0 

In the small fiow-time limit t — )■ 0, — ^ -^ +oo and the running cou- 
pling g{—C,) (I4.17P becomes very small thanks to the asymptotic freedom. 
Thus, the right-hand sides of Eqs. (I4.1ip - (l4.14p allow us to compute the 
small fiow-time behavior of the coefficients by using the perturbation theory. 

4-2. Lowest order approximation and the asymptotic formula 



By substituting the solution of the fiow equation (13. ip (see Ref. J33|) in 
the tree-level approximation to Eq. (13. 7p . we have 

CT{t) = gl (4.20) 

simply because our energy-momentum tensor (12. 4p is proportional to 1/gQ- 
If we apply the right-hand side of Eq. (14.111) to this expression by substi- 
tuting Eq. (I4.17p . however, it depends on ^/8tojJo while the left-hand side 



11 



of Eq. (14.111) does not. This shows that Cxit) should depend on g'^ and y/Sifx 
through a particular combination as (for D = 4) 

CT{t) =9^{l + 2bog^ [\n{^/8i^i) + Ci] + Oig^)} , (4.21) 

where Ci is a constant. Similarly, since the lowest-order approximation 
in Eqs. (EZD and (ESD yields 



CE{t) = 1, Cs{t) = --^QoH ^' 



(4.22) 



where the latter follows from Eq. O3.10p . from Eq. (I4.14p we have 

^ = -|/ {l + 2bog' [ln(v^/x) + c,] + 0{g')] , (4.23) 

where C2 is another constantly 

Applying Eqs. (14. lip and (I4.14p to above expressions and using Eq. (I4.17p . 
we finally have the asymptotic behaviors of the coefficients in Eq. (13. lip . 



and 



and hence 



^ '^r2^-2bo 



exit) 



ln(V8tA) + ci 



csjt) 

CEif) 

csjt) 

CT{t)cE{t) 






2 -26o [ln(y8tA) + ca] 



t^o+ 



&0 



Cl -C2 



-ln(V8tA) 



(4.26) 

(4.27) 
(4.28) 



^Using Eq. (|T^ in Eq. ((XTU| . we have 



C5(0 = -y5Nl + 2W 



ln(y8i/x)+ci + -^ 



Oig' 



and then using Eq. (|4.23l) . 



CB(t) = l + 26o5^ C1-C2 + 



26g 



0{9% 



(4.24) 



(4.25) 
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That is, 



6^,[E{t,x)-{E{t,x))] 



1 _ Ci - C2 



2 L -ln(V8tA). 

(4.29) 

This is the relation that we were seeking: One can obtain the correctly- 
normahzed conserved energy-momentum tensor from the small flow-time be- 
havior of gauge invariant products given by the Yang-Mills gradient flow. 
It is interesting to note that the leading t — )■ behavior is completely in- 
dependent of the detailed definition of the gradient flow; the structure and 
coefficients follow solely from the finiteness of the local products and the 
renormalizability of the Yang-Mills theory. The sub-leading corrections in 
the asymptotic form, i.e., the coefficients ci and C2, depend on the detailed 
definition of the gradient fiow; in the next section, we compute the constants 
Ci and C2 and we have 

7 

ci = In y^ + — ~ 1.00986, (4.30) 

16 

C2 = lnv^+^-^ + ||^ 0.812034. (4.31) 



Finally, a possible method to determine the factor ln(v8tA) in Eq. (I4.29p . 
i.e., the fiow time t in the unit of the one- loop A parameter (14.181) . for small 
fiow times is to use the expectation value of the density operator, Eq. (13. 6p . 
For this quantity, by applying Eqs. (14. 13p and (I4.17P to the result of one- 
loop calculation, Eqs. (2.28) and (2.29) of Ref. |3l| (specialized to the pure 
Yang-Mills theory), we have the asymptotic form. 



.2/z:.,. ,, t^o+3(Ar2 



where 



t' (E(t, x)) ~" ^^^ ^ -j= -, (4.32) 

^ ^ ^^ 1287r2 -2bo [ln(y8tA) + c] ' ^ ^ 

26 9 

c = ln(2v^) + In 3 -1.60396. (4.33) 



One may use this asymptotic representation for ln(A/8tA) in Eq. (14.29 



^Practically, one will use Eq. (|4.29|) to compute i^{r^,y}i?(a;) from t^Uf^,y{t,x) 
and t^E{t, x). Then, from the value of V&tA, one can deduce {Tfj,„} n.{x) / A'^ . 

13 



5. One-loop calculation of coefficient functions 

For calculational convenience, we define coefficient functions F{t) and G{t) 



by 



G';^it,x)G:^it,x) - {G;^it,x)G:^it,x)) 



= F{t) {F^,F^%}^ (x) + G{t)6,. {f;^F;,}^ (x) + 0{t). (5.1) 

Eq. (13. 5p then becomes (for D = 4), 



U^,{t,x)=F{t) 
= F{t) 



rpa pa I 
L tip up} 



MP i^PJ R 



[X 



^S^JpX {Fl^aFxaJn 



+ 0{t) 



1 



{f::pF:,}^ (x) - -6,, ^1 - ^ ) {f;,f;J^ 



(3 



x] 



+ 0(t), 
(5.2) 



wliere we liave used Eq. (I2.19p . Rewriting this in favor of the energy- 
momentum tensor (12.121) with Eqs. (12.131) and (I2.14p . we have 



U^,{t,x) = F{t) 
= F{t) 
Comparison with Eq. (13. 7p then shows 



{f^ppF:,}^{x)--6,^6,,{f;^^fi} 



R 



[X 



+ 0(t) 



/{T,.}«(x) + ^v{^;.^;.} 



/3 



R 



[X 



+ 0{t). (5.3) 



CT{t)=g'F{t), cs{t) = ^F{t). 



•29 



(5.4) 



Similarly, for Eq. (13. 8p . 



E{t, x) = {E{t, x)) + -F{t)6,, {f;^fi}^ (x) + G{t) {f;,f;J^ (x) + 0(t) 



4 



/3 



1 - ^ ) F{t) + 4G(t) 



_pa pa 
4 po- per 



(x) + 0(t). 



R 



and therefore 



CE{t) 



(3 



l--]Fit)+AGit). 



(5.5) 



(5.6) 
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This implies, for the ratio fl4.23p . 

cs{t) (3 1 



CEit) 



2gl-l + AG{t)/F{t) 



23 



4 + 26§ 



AG{t) + 0{g' 



(5.7) 



To find the coefficient functions F{t) and G{t) in Eq. (15. ip . we consider 
the correlation function 



{Gl^{t,x)Gl^{t,x)A\{w)A^^{v)). 



(5J 



To 0(5'q), there are 17 ffow-line Feynman diagrams (Figs. [T HT7I) that con- 
tribute to this correlation function. In the figures, gauge potentials at the 
flow time t, B^{t,x), are represented by little filled squares; the open circle 
denotes the fiow-time vertex and the full circle is the conventional vertex in 
the Yang-Mills theory. We refer the reader to Ref. 33| for the details of the 
Feynman rules for the fiow-line diagram. 

To read off the coefficient functions F{t) and G{t) in Eq. (15. ip from the 



correlation function (15. Sp . we consider the vertex functions, that is, ampu- 
tated diagrams in which the external propagators of the original Yang-Mills 
theory are truncated. Therefore, the diagrams Figs. [101 [12] and [171 those pro- 
vide only the conventional wave function renormalization, should be omitted 
in the computation of F{t) and (^(t)^ On the other hand, the fiow-line 
propagators |33|, the arrowed straight lines in the diagrams, should not be 
truncated because these are not propagators in the quantum field theory but 
instead represent the time evolution along the fiow time. 
The tree-level contribution to the vertex function is, 



Fig.[I]=5 



pa 






Fl{x)Fl{x) + 0{t), 



(5.9) 



®More precisely, those diagrams are different from the conventional Feynman diagrams 
in that the propagators carry an additional factor e"**' (in the Feynman gauge), where p 
is the external momentum. This factor is however irrelevant in the present computation 
of the coefficients of operators with lowest number of derivatives. 
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Figure 9: 
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Figure 16: 



Figure 17: 
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where 



M^) 



e'^^i^lp), 






and, here and in what follows, the alternating-sign symbol implies. 



tfj,pua i (/i ^ P, Z^ ^ O") = t^pjyo- — t 



pfiua '^fipav ~r <Jppav 



(5.10) 



(5.11) 



This tree- level result had been used in obtaining Eq. (I4.20p . 
The vacuum expectation value in the lowest order is 



3 



-T:^e 



-2tP , 



^p^V^ pa 



Sp^ ±{^1^ p,u ^ a) 



128712^^' "^)^o^'^^^- 



1 



(5.12) 



Now as an example of the computation of one-loop flow-line Feynman 
diagrams, we briefly illustrate the computation of Fig [131 A straightforward 
application of the Feynman rules in Ref. 33| in the "Feynman gauge" in 
which the gauge parameters are taken as Aq = ao = 1 yields the expression. 



Fig.m=Ngl6, 



0"pcr 



Ap+'i)-A''^{p)A''Jq) 



p,q 

1 1 1 



. Q-t(p+i)\~t{q-l)^ 



, {p + ey p {q - if 

X i{P + ^)m^(? - ^)u 

X [6px{-p - 2£)« + 6xa{i - p)p + 6c,p{2p + £)a] 

X [5a.(-2£ + q)^ + 5„p{-2q + £)a + 5px{q + i)a] 



(5.13) 



To find the coefficients F{t) and G{t) in Eq. (15. ip . we write this vertex 
function as 

' e'(^^'^>Al{p)Al{q)Mp,,^p{p,q), (5.14) 



P.9 



and find the coefficients of 
and 



-2p-q5^Jai3, 



(5.15) 
(5.16) 
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respectively in M^y^apiVi l)- Fo^' this, we first exponentiate the denominators 
in Eq. fl5.13p by using 



OO /"OO 



(p + ^Y {q -^y Jo Jo 

We then simply expand the integrand with respect to the external momenta 
p and q to 0{p,q). The fiow-time evolution factor e~^*^ in the integrand 
makes the integral f l5.13p UV finite for any dimension D. On the other hand, 
there always exists a complex domain of D such that the integral is infrared 
finite; this provides the analytic continuation of the integral such that 

P (47r)«/2 D/2 - 1 ' ^ ^ 



-^/'V> (5-19) 



J_ -a^2 _ 1 1 

1 _.2„ „ „ „ 1 i_ 

(47r)^/2 2(D + 2) 

(5.20) 



■e'^'^'^iiii = = = a^^^^'^(6 6 +6 6 +6 6) 



— e^"^^£ f £ £ £ i 



a "^'^ ^ ('^/xi/'^po-'^a/s + 14 permutations) . (5-21) 



(47r)^/2 4(1} + 4) 

Then it is straightforward to find the coefficients of Eqs. (I5.15P and (I5.16P 
which directly give the contribution to the functions F(t) and Gii). 

In Table [H we summarize the contribution of each diagram computed in 
the above method in the unit of 

^«. (6.22) 

In the last line of the table, "Z factors" implies the contributions of the one- 
loop operator renormalization factors, Zt (I2.13P and Zm (I2.22p through the 
tree- level diagram, Eq. (15. 9p (recall Eq. (I2.10p ). We see that those operator 
renormalization factors precisely cancel the residues of 1/e and make the 
coefficients F(t) and Gii) finite; this is precisely what we expect from the 
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Table 1: The contributions of each flow-hne Feynman diagram (in the Feynnian gauge) 
to the coefficients of Eqs. (|5.15|) and (|5.16p respectively in the unit of Eq. (|5.22[) . These 
correspond to the coefficient functions F{t) and G{t) in Eq. ()5.ip . 



F{t) G{t) 



Fig. [3] 
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general argument. From the results in the table, we then have 

(5.23) 



F{t) = 1 + 2bog^ 
G{t) = -\h,g' 



ln(y8t^)+lnv^+ — 
lo 



ln(y8t/i) + lnv^+ — 



(5.24) 



Finally, comparison with the formulas (15. 4p . (15. 7p . (I4.2ip and (I4.23P shows 
the results quoted in the previous section, Eqs. (I4.30p and (14.3 ip . Note that 
the coefficients of ln(-\/8t/i) in the explicit one-loop calculation (Eqs. (I5.23P 
and (I5.24P ) are in agreement with those by the general argument on the basis 
of the renormalization group equations and the trace anomaly (Eqs. (I4.2ip 
and (14.231) ). This agreement provides a consistency check of our one- loop 
calculation and supports the correctness of our reasoning. 

6. Conclusion 

In the present paper, we derived a formula that relates the short flow-time 
behavior of some gauge invariant local products generated by the Yang-Mills 
gradient flow and the correctly-normalized conserved energy-momentum ten- 
sor in the Yang-Mills theory. Our main result is Eq. (I4.29p . The right-hand 
side of Eq. (I4.29p can be computed by the Wilson flow in lattice gauge the- 
ory with appropriate discretizations of operators, Eqs. (13. 5p and (13. 6 p (see 



Refs. [31|, |3J], for example). Here, the continuum limit a — )■ must be 
taken ffist and then the t — ;■ limit is taken afterwards; otherwise our basic 
reasoning does not hold. 

Although the formula (I4.29P should be mathematically correct, the prac- 
tical usefulness of Eq. (I4.29P is a separate issue and has to be carefully ex- 
amined numericallylj Since the lattice spacing a must be sufficiently smaller 
than the square-root of the flow time \/8t for our reasoning to work, the 
reliable application of Eq. (I4.29P will require rather small lattice spacings. 
One also worries about the contamination by higher dimensional operators 
(i.e., 0{t) terms in Eqs. (13. 7p and (13.81) ) and the finite-size effect which we 
did not take into account in the present paper. If our strategy turns to be 
practically feasible, it provides a completely new method to compute correla- 
tion functions containing a well-defined energy-momentum tensor. It is clear 



^We hope to return this problem in the near future. 
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that the present approach to the energy-momentuin tensor on the lattice is 
not hniited to the pure Yang-Mills theory although the treatment might be 
slightly more complicated with the presence of other fields. The application 
will then include, the determination of the shear and bulk viscosities (see 



Refs. |43|, 44, for example), the measurement of thermodynamical quantities 



(see Ref. [45| and references cited therein), the mass and the decay constant 
of the pseudo Nambu-Goldstone boson associated with the (approximate) 



dilatation invariance (see Ref. [46| and references cited therein), and so on. 
It is also clear that our basic idea that operators defined with the lattice 
regularization and in the continuum theory can be related through the gradi- 
ent fiow is not limited to the energy-momentum tensor. It might be possible 
to construct, for example, an ideal chiral current or an ideal supercurrent 
on the lattice, from the small fiow-time limit of local products. It appears 
interesting to pursue this idea. 
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